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Abstract 

Bonnor stars are regular static compact configurations in equilibrium, composed of an extremal dust fluid, 
i.e., a charged dust fluid where the mass density is equal to the charge density in appropriate units and up to 
a sign, joined to a suitable exterior vacuum solution, both within Newtonian gravity and general relativity. 
In four dimensions, these configurations obey the Majumdar-Papapetrou system of equations, in one case, 
the system is a particular setup of Newtonian gravity coupled to Coulomb electricity and electrically charged 
matter or fluid, in the other case, the system is a particular setup of general relativity coupled to Maxwell 
electromagnetism and electrically charged matter or fluid, where the corresponding gravitational potential 
is a specially simple function of the electric potential field and the fluid, when there is one, is made of 
extremal dust. Since the Majumdar-Papapetrou system can be generalized to d spacetime dimensions, 
as has been previously done, and higher dimensional scenarios can be important in gravitational physics, 
it is natural to study this type of Bonnor solutions in higher dimensions, d > 4. As a preparation, we 
analyze Newton-Coulomb theory with an electrically charged fluid in a Majumdar-Papapetrou context, in 
d = n + 1 spacetime dimensions, with n being the number of spatial dimensions. We show that within the 
Newtonian theory, in vacuum, the Majumdar-Papapetrou relation for the gravitational potential in terms of 
the electric potential, and its related Weyl relation, are equivalent, in contrast with general relativity where 
they are distinct. We study a class of spherically symmetric Bonnor stars within this theory. Under sufficient 
compactification they form point mass charged Newtonian singularities. We then study the analogue type 
systems in the Einstein-Maxwell theory with an electrically charged fluid. Drawing on our previous work 
on the d-dimensional Majumdar-Papapetrou system, we restate some properties of this system. We obtain 
spherically symmetric Bonnor star solutions in d = n + 1 spacetime dimensions. We show that these 
stars, under sufficient compactification, form d-dimcnsional quasi-black holes. We also show that in the 
appropriate low gravity limit theses solutions turn into the solutions of Newtonian gravity, i.e., they are 
quasi-Newtonian Bonnor stars. In this connection, we note that the star solutions in Majumdar-Papapetrou 
Newtonian gravity, when contrasted to those solutions in Majumdar-Papapetrou general relativity, display 
clearly the branching off of the high density objects that may arise in the strong field regime of each theory, 
mild singularities in one theory, quasi-black holes in the other. Another important feature worth of mention 
is that, whereas there are no solutions for Newtonian or relativistic stars supported by degenerate pressure 
in higher dimensions, higher dimensional Bonnor stars, supported by electric repulsion, do indeed have 
solutions within Newtonian gravity and general relativity. So the existence of stars in higher dimensions 
depends on the number of dimensions itself, and on the underlying field content of those stars. 

PACS numbers: 04.50Gh, 04.40Nr, ll.lOKk 
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I. INTRODUCTION 



A. Definition 



Extremal charged dust, or simply extremal dust, is understood as charged dust fluid, or matter, 
with the mass density being equal to the charge density, in appropriate units, which implies that for 
each such a dust particle, eventually composing a system, the gravitational attraction is precisely 
balanced by the electric repulsion, both within Newtonian gravity coupled to Coulomb electricity 
and an electrically charged fluid or matter, i.e., the Newton-Coulomb with charged fluid system, 
and within general relativity coupled to Maxwell electromagnetism and an electrically charged 
fluid, i.e., the Einstein-Maxwell with charged fluid system. Bonnor stars are then defined as 
regular static equilibrium configurations, in Newtonian and general relativity contexts, composed 
of extremal dust, with a finite boundary appropriately attached to an asymptotically flat regular 
extremal charged vacuum, and where the configuration of the matter dust can have any shape, a 
spherical symmetric shape being of special interest, due to the added symmetry and due to the 
fact that it can be joined to an asymptotically flat regular extremal outer Reissner-Nordstrom 
spacetime in the general relativistic case, with mass M equal to charge Q in appropriate units. 
Bonnor stars appear in d = n + 1 spacetime dimensions, where n is the number of the spatial 
dimensions. The initial studies where performed for d = 4. Bonnor stars have also been called 
Majumdar-Papapetrou stars, but here we reserve the name Majumdar-Papapetrou for the type 
of matter, and name the whole system, namely Majumdar-Papapetrou matter plus vacuum plus 
junction, Bonnor star. 



B. Four dimensional analyses 



1. Context 



Such stars were studied mainl y w ithin general relativity, although with some incursions onto 
Newtonian gravity, by Bonnor and in several other works, see e.g., [ul [bJ] for Bonnor stars 

properly said, and in [15[] for a variant where there is no need for a junction. One striking property 
of these solutions found in [ij], Q, [ijj is that when they approach their gravitational radius in a 
static sequence of configurations, these stars do not form black holes, but rather quasi-black holes, 
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where a quasi-black hole is an object indistinguishable to the exterior from a black hole but with 



different intrinsic properties. In 



16 



17 



the properties of quasi-black holes formed from Bonnor 



stars were studied in detail. See also [18(] for a further study on the properties of Bonnor stars. 



2. Vacuum Majumdar-Papapetrou solutions 

Within Newtonian gravity coupled to Coulomb electricity, i.e., the Newton-Coulomb system, 
electrically charged solutions in vacuum represent charged point masses. Within general relativ- 
ity, electrically charged solutions, in vacuum, have to be analyzed through the Einstein-Maxwell 
system of equations, where one couples Einstein gravity to Maxwell electromagnetism. These so- 
lutions were found just after general relativity was formulated. On one hand, Reissner [ijj], then 



Xonistrom [20I ] . then Jeffery [2l[ |. hit on the static vacuum charged spherically symmetric solution, 
the Reissner-Nordstrom solution, with its two parameters, the mass M and the charge Q. We now 
know, that for y/GM < eQ (we put the speed of light c = 1 throughout), where G is Newton's 
gravitational constant in four spacetime dimensions (G4 = G) and e = 1 if the charge is positive 
and e = — 1 if the charge is negative, one has a charged naked singularity, for yGM > eQone has 



221 for an 



a Reissner-Nordstrom black hole, and V GM = eQ one has an extremal black hole (see 
early discussion, and [23] for a full discussion of this type of solutions). 

On the other hand, with the purpose of seeking vacuum static solutions electrically charged, a 



different route was originated from Weyl [24j |. the route that will take us to the Bonnor stars [ll|-[l8]. 
Wanting to go a step further from spherical symmetry he sought axial symmetry. Define the metric 
component goo as 500 — W 2 , where, depending on the situation, it can be more convenient to define 
U = W~ l , i.e., goo = U~ 2 . Then, Weyl asked himself, within Einstein-Maxwell theory, what would 
happen if W 2 , in a static axisymmetric vacuum electric system, is to have a functional dependence 
on the electric potential field ip, i.e., the Weyl ansatz goo = dooip) or equivalently W = W(ip). 
He found first what is now called the Weyl relation, i.e., W 2 = (ao — ey/Gip^j + &o> where ao 
and 6q are constants of integration and G is Newton's gravitational constant, and second that the 



25] 



spatial components of the metric had to obey other specific differential equations. Majumdar 
made several improvements on Weyl's work. He showed that the Weyl relation, if it existed was 
independent of the symmetry, axial or otherwise. But further, he showed, still in vacuum, that 
if the relation was to be a perfect square, so that, W = ao — eyGtp, then the spatial part of 
the metric could be put in a simple form, as 1/W 2 , i.e. U 2 , times the flat spatial metric, and 
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the Einstein-Maxwell system of equations would reduce to one single equation for W, i.e. for 
U, a Laplace equation in flat space. In this perfect square case, one can show that specializing 
to spherical symmetry, the mass M of the solution is equal to its charge Q, yGM = eQ. This 
makes contact with the Reissner-Nordstrom family of solutions through the extremal solution, 
\[GM = eQ, although not through the other ones, since the Reissner-Nordstrom family generically 
does not admit a functional relationship between the metric and the electric potentials. These 
vacuum yGM = eQ solutions were further commented by Papapetrou 26], who also understood 



that since the gravitational attraction is equal to the electric repulsion for such objects one could 
have many discrete such objects scattered at will in space, that it would also give a vacuum 
static configuration solution, with no symmetry whatsoever. The perfect square relation is usually 
called Majumdar-Papapetrou relation, as we do here, although sometimes it is called, perhaps 
more appropriately Weyl-Majumdar relation. The complete understanding of a single extremal 
Reissner-Nordstrom, also a Majumdar-Papapetrou solution, solution was achieved by Carter j^ . 
through a Carter-Penrose diagram, and the complete understanding of the vacuum Majumdar- 
Papapetrou solutions, with many extremal black holes scattered around was performed by Hartle 
and Hawking [2^, who have done the maximal analytical continuation in the molds of Carter [27]. 



3. Beyond vacuum: dust Majumdar-Papapetrou solutions 

Things become more interesting if one goes beyond vacuum and puts matter into the Newton- 
Coulomb system of equations and into the Einstein-Maxwell system of equations. Majumdar 



251 ] and Papapetrou 26] understood this, and within general relativity, showed, for some special 
restrictions on the metric inspired from the vacuum case, such as the relation W = clq — e\J~G p 
(which in this case can be considered an ansatz), one could find the system of equations yield a 
single equation that moreover reduces to a Poisson equation, and in which the mass density p m 
times vG is equal to the charge density p e , up to a sign, y/Gp m = ep e , with again e = dbl. That 
is, the matter is made of an extremal dust fluid. This is the Majumdar-Papapetrou condition. 
Note that the relation between the potentials we call Majumdar-Papapetrou relation, whereas the 



25, 



26(], it 



relation between the densities we call Majumdar-Papapetrou condition. As seen in 
is remarkable that a simple obvious fact in Newton-Coulomb theory with an electrically charged 
fluid, that if the mass density and charge density are equal (in appropriate (geometric) units 
where G = 1) then there is exact balancing of the gravitational and electric forces throughout the 
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matter and so there is a static solution, also holds in Einstein-Maxwell theory with an electrically 
charged fluid, with no need for further stresses, such as pressure or tension. The basic feature of 
Majumdar-Papapetrou systems is that they describe static spacetimes filled either with extremal 
charged vacuum or extremal charged dust fluids, such that the metric and electromagnetic fields 
may be characterized by two scalar functions, namely, the redshift metric function W, i.e. U , 
which plays the role of the gravitational potential, and the electric potential tp, which in turn 



obey the Majumdar-Papapetrou relation, W = on — eyGtp. Further interesting developments were 
achieved by Das [29J ]. and De and Raychaudhuri 30], who considered charged dust distributions in 
equilibrium, the way envisioned by Majumdar and Papapetrou, and showed some other conditions 
related to the functional form of the metric in terms of the electric potential and the equality 
between mass and charge densities. Das 



291 revealed that the equality between the densities 

□ 

implies the functional form on the potentials, and De and Raychaudhuri [30J proved that given 
the functional form above, and provided there are no singularities in the distribution, the equality 
of mass and charge densities follows directly from the field equations. There are other interesting 
properties of Majumdar-Papapetrou systems in the context of conformal static charged solutions 



4- Bonnor stars: junction of dust with vacuum Majumdar-Papapetrou solutions 

When one joins smoothly, within Newtonian gravity as well as within general relativity, 
Majumdar-Papapetrou interior matter solutions to Majumdar-Papapetrou exterior vacuum so- 
lutions, i.e^ the two types of solutions mentioned in the previous paragraphs, one obtains Bonnor 
stars [lj]- and their developments Bonnor stars could instead be called Majumdar- 

Papapetrou stars as was done in jl41 ]. but it is more proper to characterize the matter part as a 
Majumdar-Papapetrou system, and this combined with placing a boundary and the corresponding 
junction to a vacuum, bringing together a whole lot of new properties, as a Bonnor star. Through- 
out his works, Bonnor gradually improved the understanding of the properties of these stars. 

In the first two papers Bonnor worked out aspects of electric Majumdar-Papapetrou 

solutions in an axial symmetric vacuum and extended these results through dualities to magnetic 
solutions. In [3] a pre-Bonnor star is developed, and it is noted that for \^Gp m = ep e , in a 
Majumdar-Papapetrou system, the gravitational mass of the system is equal to the matter mass 
because the negative gravitational self-energy of the distribution is balanced by the corresponding 
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positive electrical self-energy, also pointing out that v GM = eQ models can have various interests 

n rl 

and applications. In [4, 5] Bonnor understood for the first time that although delicate, the balance 
can exist, an atom stripped off of an electron immersed in about 10 18 atoms is enough, and that the 
charge density can play an important part in the equilibrium of large bodies, further suggesting 
that it may halt gravitational collapse, at a time where large bodies studies were in vogue due 
to the appearance of quasars. It is also mentioned that bodies of arbitrary shape composed of 
such extremal dust can exist, constructing explicitly a spherically symmetric solution, the first 
Bonnor star. The way it is constructed delineates a standard way of finding such type of solutions. 
Assuming a given form for the gravitational potential U one can find the density distribution, 
and one hopes that that assumption yields a physical distribution of charged dust matter. It is 
not a method for solving the differential equation of the Majumdar-Papapetrou problem, it is an 
art of correct guessing. In fl fl, both w„r ts in collaboration with Wich™ri y a, interesting 
physical properties of some Bonnor stars are discussed. In particular in [g] the name electrically 
counterpoised dust, or ECD for short, is coined for the first time for \[Gp m = ep e dust, i.e., for 
what we call and will always call, less contrively perhaps, extremal charged dust, with the same 
acronym. Spherically symmetric exact solutions are studied with the virtue that even when the 
solutions are about to form a horizon the energy density p m is finite. Prolate solutions are also 
studied showing that in the disk solution limit the energy density p m is infinite, naturally. Also, 
in particular, in [7| several important attributes of the solutions are perceived. First, it is noticed 
that, although no doubt, matter thus delicately balanced is rare, it is physically possible and 
easily understood. Second, it is shown that one can construct spheres of matter where infinite 
redshifts of light emanating from the surface are attainable, whereas in an interior Schwarzschild 
solution, say, only finite redshifts are possible. Third, it displays exact solutions for spheroidal 
configurations, and mentions that near the spheroidal horizon, non-spherically symmetric features 
are filtered out. Fourth, it is argued convincingly that these solutions are stable. In [§] the study 
is interesting with queer results. First, there is an incursion into solutions of the Newton-Coulomb 
with an electrically charged fluid theory, where it is shown that for a given Newtonian potential, 
call it V, there are equilibrium non Bonnor star solutions, not obeying the Majumdar-Papapetrou 
matter condition, although these are singular. It also shows the analogue of De and Raychaudhuri's 
theorem [3o| . i.e., that Newtonian systems which do not obey the Majumdar-Papapetrou condition, 
of the equality of mass and charge densities, and which have equipotential surfaces, are singular. 
Second, in turning into general relativity, with axial symmetry, spacetimes obeying the Majumdar- 
Papapetrou condition are found, one of them being of physical interest with positive energy density 
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p m , and the others of less interest. In [{], [l(| it is understood, perhaps for the first time, that 
when the radius of the configuration r$ approaches the horizon radius, i.e., ro = M, where r is 
the Schwarzschildean radial coordinate and M the mass of the configuration, the spacetime is 
somehow singular, being thus an idea precursor of the concept of what a quasi-black hole is. In 
these works the hoop conjecture is discussed and some lower bounds in connection to it are given. 
In 



111 ], spheroidal bodies made of extremal charged dust are studied in connection still with the 
hoop conjecture and also with the isoperimetric conjecture, which says that under certain general 
conditions M > {A/\Q-k) 1 I 2 , where A is the area of the apparent horizon, and M the mass of the 
configuration. In [12|] it is reinforced that spherically symmetric configurations made of physically 
reasonable matter, though admittedly not widely available, i.e., made of extremal charged dust, 
yield solutions that can come arbitrarily close to the horizon of an extremal black hole, and a 
general class of such solutions is constructed by correct guessing. 



Bonnor stars were studied further by other authors. In [13] a thick shell solution of Bonnor type 
was found. In []jj] it was noted that Bonnor star solutions and gravitational magnetic monopole 
solutions have striking similar properties, and a comparison of both solutions was performed and 
discussed thoroughly. Previously, Lemos and Weinberg [l5l |. seeing in Bonnor 's wake that these 
stars, which are made of normal matter obeying the several important energy conditions, can 
probe deeply the spacetime structure, proposed new solutions, extended Bonnor star systems with 
a more sophisticated density distribution asymptotic to an extreme Reissner-Nordstrom solution, 
not needing any junction. Similar properties were found for Bonnor stars properly said as well as 
for extended Bonnor stars. Most notably, is the fact that at the threshold of the formation of an 
event horizon the system displays a very peculiar trait, instead of an extremal black hole one has a 
quasi-black hole, with the formation of a quasihorizon instead of the usual event horizon. Although 
to external observers the system looks like an extremal black hole, its internal properties are very 
different from what one could expect in the case of a standard black hole. These properties, along 
similar ones of gravitational magnetic monopoles and glued vacuum solutions with shells, have 
been analyzed in [it]. \\\. In 18] other attributes of these systems were explored. 



5. Further connections 



One can associate these Bonnor stars to several related topics, (i) Both, astrophysically and 
physically, Bonnor stars are of interest. On one hand, they can be realized if a gravitating sphere, 
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of neutral hydrogen which has lost a fraction 1CP 18 of its electrons, forms. On the other hand, 
they are supersymmetric solutions of N = 2 supergravity [32|, so are of interest in an elementary 
particle context, (ii) A matter always of maximal interest is the stability of the systems one is 
considering, in this case, the Bonnor stars. Interestingly enough it was found, through different 
methods, that these stars are neutrally stable. Firstly, Omote and Sato [33J found this stability 
criterion using both an energy method and a small adiabatic radial oscillation method, results 



which were later confirmed in 



34 



351 ] . (iii) When discussing static equilibrium configurations it is 
always important to discuss the Buchdahl limits, where for instance for a perfect fluid sphere, the 
star cannot reach beyond ro < 9/8rs c h w > where tq is the star radius, rs c hw is the Schwarzschild 
radius, rs c h w = 1GM [36|], and r is the Schwarzschildean radial coordinate. Contrarily, for Bonnor 
stars, stars made of extremal charged matter, the limits are precisely the horizon radius as was 
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40 



first noted by Bonnor , and then in subsequent works 

for interesting discussions on the Buchdahl limits for charged stars, (iv) The hoop conjecture is 



13j-[17t], see also 



39, 



4l| 



charge 



relevant for these stars as was first noticed by Bonnor [9j, [l(| (see also [111]). The conjecture states 
that a black hole forms when matter of mass M is compacted within a given definite hoop, in 



421 ] taken to be ~ 4ttGM. Later, it was shown that the hoop should be reduced for extremal 
charged matter to ~ 2irGM [9j, lQ l. However, it seems that systems like Bonnor stars violate 

n n 

it, since no black hole forms ever, only a quasi-black hole [la . I X 71 ] - (v) A pertinent question, 
specially related to stars, is whether they can form from gravitational collapse or not. The issue 
of the collapse of extremal charged dust solutions has not been studied in detail, see, however, the 
interesting work of De 4s| . (vi) Concerning the generalization of Bonnor stars to include pressure 
terms, and thus go beyond dust matter there are some stimulating developments. For instance, 
still within the Majumdar-Papapetrou ansatz for the potential, W = ao — e^/~G ip, systems with 



pressure were studied by Ida [44|, where one finds, with an additional ansatz for the pressure, a 
Helmholtz type equation which can be solved, in the case the pressure is zero see also These are 
thus extensions of Bonnor stars, stars that include charged matter, non-extremal, and pressure. 
Extensions to solutions with potentials different from the Majumdar-Papapetrou potential, and 
even different to Weyl's potential, were done in These solutions include pressure and 

have interesting structure. Charged stars with pressure, were studied numerically in [48], a paper 
that has attracted some attention, where the limiting configuration is found to have mass equal 
to charge, in appropriate units, being thus a Bonnor star. In [49] a set of static charged solutions 
with pressure were studied and it was proposed that their gravitational collapse would lead to the 
formation of a charged Reissner-Nordstrom black hole, (vii) There are many other solutions of 
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charged matter in various situations which have been discovered throughout the years. Many of 
them are interesting and would deserve a review, but there are too many to be quoted here, see 
for a very partial list, (viii) Charged gravitating solutions have been also used to study Abraham 
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52 



53|, and also 



,3, 



type models for the electron, with and without Poincare stresses see, e.g, 
and [17]. (ix) A related issue to Bonnor star solutions and quasi-black holes is the set of black 
holes devised by Bardeen [5J], in which the interior to the horizon is nonsingular. These solutions 



are magnetically charged, instead of electrically charged, and have been further explored in 
The connection with the quasi-black holes is that there is a theorem by Borde [56] which says that 
if there is no singularity inside the event horizon then the regular solutions have different inside 
and outside topologies. Now, it is not possible to put extremal charged dust, with positive rest 
mass, inside an extremal black hole, a la Bardeen, a result first found in <i-dimensional studies (see 
below). So physical (positive rest mass) Bonnor stars do not provide Bardeen like solutions. On 
the other hand quasi-black holes are not true black holes, but have a weird topology and properties 
[l^.[l7]. approaching considerably the topology change of Borde. For further connections of Bonnor 
stars and quasi-black holes with other issues, such as no hair theorems, naked black holes, objects 
that mimic black holes, and the entropy issue, see [nl [l7|. 



C. Higher dimensional analyses 



1. Context 



The possibility of the existence of extra dimensions arise in several theoretical schemes. In what 
is called a Kaluza-Klein unification model, the unification idea has emerged first as a way of unifying 
the gravitational and electromagnetic fields in five spacetime dimensions, and later the gravitational 
and Yang-Mills fields in seven spacetime dimensions. Within this idea the gravitational field 
in higher dimensions gives rise to the gravitational field itself and the other possible fields in 
four dimensions. Later, the Kaluza-Klein process was enforced in theories which start from the 
outset in higher dimensions, such as supergravity or string-M theory, which can have up to eleven 
dimensions. In the course of reducing the dimensions to four, a profusion of new fields materialize, 
see {57] ] for the original papers. These schemes require that the extra dimensions are compactified 
in Planck size manifolds, and so due to the lack of a properly accepted theory at these scales 
it is very hard to do physics on the extra dimensions. It has now appeared an idea that makes 
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the higher dimensions large, when compared to the Planck scale, which means, if correct, it can 
have measurable consequences on current or near future experiments. By postulating that the 
gravitational field propagates also in at least extra three space dimensions, while electromagnetism 
and the standard model fields propagate only in our universe, the brane, it is possible to reduce 
the Planck scale to the electroweak scale and make the extra dimensions large, of the order of 
hundredths of a centimeter or a bit less. The hierarchy problem, of understanding the huge 
differences in the gravity and electroweak scales, is now pushed into the acceptance of the large 



extra dimensions (see 



58, 



59 



60|, see also [6l( for possible developments) 



Now, since within this arrangement gravity is an electroweak scale phenomenon, so are black 
holes. Thus, for instance, by splashing electrically charged particles together black holes or other 
gravitational objects can be created in higher dimensions with the charge remaining in the brane. 
In scenarios with extra dimensions it is thus important to study charged solutions in connection 
to the formation of these tiny black holes because the charge and the solutions suggest that the 
charge may halt gravitational collapse. Solutions for charged objects in such a frame are certainly 
not spherically symmetric, so not Reissner-Nordstrom, and at the moment have not been found. 
Nonetheless, it is certainly of interest to consider spherically symmetric electrically charged solu- 
tions in higher dimensions because, first such a study can give an idea of how the existence of 
the charge influences the solution, and second, other charged fields, analogous in many respects 
to the electromagnetic field, may propagate in the higher dimensions, making Maxwell electrically 
charged solutions prototype solutions. 

In addition, related to studies on the role played by the dimensionality of space on the laws of 
physics and its connection to the anthropic principle, it has been shown that there are no Newtonian 
solutions for stars supported by degenerate pressure in higher dimensions, i.e., a higher dimensional 
self-gravitating Fermi gas either collapses into a black hole or evaporates. Indeed, interesting papers 
discussing degenerate stars, like white dwarfs and neutron stars, in higher dimensions have appeared 



62, 



63j. In 



62[ | a complete heuristic study, following the original work of Landau (see, e.g, 



was performed. Then in 



631 ] . the full study, following the original works of Chandrasekhar (see, e.g, 
651] ). was completed. The main conclusion is that there are no Newtonian solutions for degenerate 
stars in higher dimensions, thus no general relativistic solutions also, the Fermi pressure energy 
cannot balance the gravitational energy. Of course this may not follow for other stars. Stars 
supported by classical gas pressure may perhaps exist in higher dimensions, no conclusive study 
has been presented so far. Thus, it is of interest to show whether stars, supported by electric 
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repulsion, such as Bonnor stars, do have solutions within Newtonian gravity and general relativity. 
In case there are solutions, one shows by example, that the existence of stars in higher dimensions 
depends both on the number of dimensions itself, and on the underlying field content of the stars 
themselves. 

It is thus important, for the reasons just raised, to study Bonnor stars in higher dimensions, 
prior to compactification of any sort. 



2. Vacuum Majumdar-Papapetrou solutions 

Electrically charged solutions in vacuum in d dimensions within Newton-Coulomb theory are 
a direct generalization from four dimensions. Within Einstein-Maxwell theory the d-dimensional 
solutions were found by Tangherlini 



66J , with a prescient discussion on the physical laws and their 
relationship to the three dimensionality of space. These solutions generalize the four dimensional 
Reissner-Nordstrom solutions, and they also have the mass M and the charge Q, as the higher 
dimensional parameters, such that for \/G~d M < eQ one has a charged naked singularity, for 
\/Gd M > eQ one has a Reissner-Nordstrom black hole, and y/Gj M = eQ one has an extremal 
black hole, where e = ±1 depending on the sign of the charge. Here Gd is the d-dimensional 
Newton's gravitational constant, where in four spacetime dimensions we put G4 = G (see Appendix 
[A] for more on this). If one takes Weyl's and Majumdar's route into higher dimensional Einstein- 



Maxwell theory, see now 67], and seeks the initial ansatz that the metric potential W, or its inverse 
U = W^ 1 , depends on the electric potential tp, i.e., W(ip) with goo = W 2 , one finds the relation 
W 2 = (ao — tVGd 93) 2 + bo, also independent of the symmetry. In the perfect square Majumdar- 



Papapetrou case, i.e. W = ao — ey/Gd </>, one can also show that specializing to spherical symmetry, 
the mass M of the solution is equal in appropriate units to its charge Q, x/Gj M = eQ. This makes 
contact with the Tangherlini black holes through the extremal solution \/Gj M = eQ, although 
not through the other ones, since the Tangherlini family generically does not admit a functional 
relationship between the metric and the electric potentials. The complete understanding of a single 
extremal Reissner-Nordstrom solution can also be achieved through Carter-Penrose diagrams, and 
the complete understanding of the vacuum Majumdar-Papapetrou solutions, with many extremal 
black holes scattered around in d dimensions was performed in 68]. 



Moreover, it is interesting to note that if instead of working in Einstein-Maxwell theory, one 
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works in string-M theory or in supergravity theory in eleven dimensions, there are Majumdar- 



Papapetrou type solutions, in the sense that the attraction due to the gravitational fie 



69l | , as well as [70(] , for reviews 



d is counter- 



balanced by the repulsion of the charged field of the theory, see, e.g, 

on this topic (see also [ti| for a review on black hole and other solutions of higher-dimensional 
vacuum general relativity and higher-dimensional supergravity theories). In eleven dimensions in 
string-M theory, there are two bosonic fields, the metric and the A3 form field which is a variant of 
the electromagnetic field with a corresponding charge, and one fermionic field. Thus the bosonic 
part is as simple as Einstein-Maxwell. One finds that for the solutions to be purely bosonic one has 
to have that the mass of the solution has to be equal to the A3 charge, in appropriate units. Note 
that this is the analogue of the extremality bound for Reissner-Nordstrom black holes. Solutions 
with mass equal to charge are called BPS (Bogomolnyi-Prasad-Sommerfield) spacetimes. The 
solutions are not point-like generically. They are brane like, and are called p-branes, or black p- 
branes, where a zero-brane is a zero dimensional object like a black hole, a one-brane is a string like 
a black string, a two-brane is a membrane and so on. Indeed, in string-M theory, where supergravity 
in eleven dimensions is a low energy theory, there are the M2-brane (a membrane, i.e., a two-brane 
electrically charged under A3), the M5-brane (a five-brane magnetically charged under A3), the 
wave solution or Aichelburg-Sexl metric, and the Kaluza-Klein monopole. All of these are BPS, the 
last two having momentum which is a form of charge. These solutions are best found and studied 
in isotropic, also called harmonic, coordinates, as is the case of Majumdar-Papapetrou solutions 
in general relativity. One can then have, for instance, many M2-branes scattered around, as one 
can have many black holes in the Majumdar-Papapetrou case, since the charged field force still 
balances the gravitational force. One can, in addition, combine the solutions with different charge 
type, for instance a M2-brane with a M5-brane, with no analogue in Majumdar-Papapetrou since 
here there is only one charge. Through careful dimensional reduction, Kaluza-Klein or otherwise, 
these solutions are also solutions of the reduced theories. Usually the branes in eleven dimensions 
are non singular and considered as solitonic objects. But when one reduces to ten dimensions, 
singularities in the solutions appear, in which case it is better to consider the branes as coupled to 
extremal dust, in the place of the singularities (see, e.g., making thus the consideration 

of extremal dust solutions in higher dimensions a subject of interest. It is also worth commenting 
that in string-M theory in eleven dimensions one can also perform some brane engineering, by 
adding together solutions of the same type of charge. It is common practice to stack an array 
of M2 electrically charged branes, for instance, and then take the continuum limit, or smear, the 
array in the correct direction, yielding a new brane with a new dimension, see, e.g., 0,0]. Of 
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course one can also do the same type of manipulation in Majumdar-Papapetrou general relativity. 
Draw an array of equally sparse extremal black holes on a line, smear them together correctly, and 
obtain a one dimensional extremal black string obeying the ci-dimensional Majumdar-Papapetrou 
equations. 



3. Beyond vacuum: dust Majumdar-Papapetrou solutions 

In d dimensions, as in four, things become more interesting if one goes beyond vacuum and puts 
an electrically charged fluid or matter into the Newtonian gravity or general relativistic systems 



of equations. Leaning on the general relativistic analysis of Majumdar 



251 ]. Lemos and Zanchin 



671 ] showed, for the special relation, or ansatz in this context, on the metric inspired from the 
vacuum case, i.e., W = ao — ty/Gd ^P-, that the whole system reduces to a single equation, a Poisson 
type equation, in which the mass density p m is equal to the charge density p e in appropriate units, 
\JGd p m = ep e . Thus, a basic feature of such a system is that, although being a system containing 
charged matter, it is described by the metric function W, the redshift function. The electric 
potential ip can then be found implicitly through the Majumdar-Papapetrou relation. It is also 
possible to generalize to d dimensions the theorem, proved in four dimensions in general relativity 
in [3^], that, provided the pressure is zero and there are no singularities in the distribution, the 



Z2(|. One 



Majumdar-Papapetrou ansatz W = ao — ey/Gd <p and condition y/Gd Pm = epe follow 
can then show that the d-dimensional Newtonian limit follows, with the four dimensional situation 
studied in Q] being a particular case. Also theorems with nonzero pressure can be render into 



d dimensions |72|], namely, that for perfect fluid solutions satisfying the Majumdar-Papapetrou 
condition the pressure is related to redshift function, as in the four dimensional case j^ . 



4- Bonnor stars: junction of dust with vacuum Majumdar-Papapetrou solutions 



In 67] it was proved that if the pressure is functionally related to the redshift function, which in 
turn obeys the Majumdar-Papapetrou relation for the potentials, then to have a surface with zero 
pressure, i.e., a star, one has to have the pressure equal to zero everywhere. This in turn means the 
star is a Bonnor star, with a d-dimensional Majumdar-Papapetrou interior and a d-dimensional 
extremal Reissner-Nordstrom exterior. This result is valid within both Newtonian gravity and 
general relativity. To discuss Bonnor star solutions in the spherically symmetric case is the aim of 
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this paper. We show that spherically symmetric Bonnor stars in d dimensions have a number of 
interesting properties. In Newtonian theory their mass and radius may be arbitrary and the object 
with the highest compression is a point electric mass, i.e., a Newtonian singularity. In general 
relativity the stars can yield very large redshifts and their exteriors can be made arbitrarily near 
to the exterior of an extremal charged black hole. Even in these extremal situations, many of their 
characteristics remain finite and non-trivial. These extremal kind of c?-dimensional systems are the 
quasi-black holes, possessing quasihorizons, already mentioned. 



5. Further connections 



As in four dimensions, in d dimensions one can try to associate Bonnor stars to several related 
topics, (i) Bonnor stars, or something related, in higher dimensions are of interest in situations 
prior to compactification. Since astrophysically, the world is already compactified to four space- 
time dimensions, the main interest in these solutions is for high energy physics, for instance in 
a large extra dimension scenario, where charged configurations in higher dimensions can be of 
interest. It would be of interest to know whether cf-dimensional Bonnor stars, for generic d, are 
supersymmetric solutions when embedded in some supergravity theory, (ii) Of course, the study 
of the stability of these higher dimensional stars is important, although we do not do it here, (iii) 
Buchdahl limits in higher dimensions have not been found neither for uncharged nor for charged 
stars. We are preparing such a study, (iv) As far as we know, there is no discussion of the hoop 
conjecture for objects in d dimensions, (v) In d dimensions it is also important to understand if 
the configurations under study can form from gravitational collapse. Collapsing and static charged 
shells in d dimensions within Einstein-Maxwell theory with an electrically charged fluid have been 

umdar-Papapetrou 



analyzed in [73f] . Static shells, with vanishing pressure, in this context are Ma. 
solutions. Collapsing charged shells in Lovelock theory have been studied in 



74 ] . (vi) One can also 



use the Maj umdar-Papapetrou relation for the potential W = oo — eyGj ip, and study systems 
with pressure in much the same way as Ida [44]. We will not discuss this type of solutions in d 
dimensions, (vii) There are a few other, non Majumdar-Papapetrou type, solutions of charged 
matter, see, e.g, the interesting ones discussed in 13]) w here charged spheres with specific 
distributions of matter, charge, and pressure were found, (viii) Electron models in the molds of 
Abraham and Poincare have not been studied in d dimensions, (ix) It would be interesting to 
study Bardeen models and Borde's theorem in d dimensions. An interesting result first derived in 
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[73| is that for a shell in d dimensions with positive proper mass there is no static solution inside 
the event horizon, the result being valid if four dimensions as mentioned above, as well as in d > 4. 
This in some sense connects with Borde's theorem 56]. 



D. Lay out 

We start by analyzing, in Section HH the Newtonian theory for charged fluids in higher dimen- 
sions, looking for static solutions. We verify in subsection III Al that if the condition y/Gj Pm = £Pe, 
with e = ±1, is to be satisfied, then there are equilibrium Bonnor star solutions in d = (n + 1)- 
dimensional spherically symmetric spacetimes, where n is the dimension of the space, see subsection 
III Bt Bonnor stars of Majumdar-Papapetrou general relativity are studied in Section IIII[ In sub- 
section IIII Al we write the basic equations and particularize them for spherical symmetry. Part 
of the subsection is devoted to review the main properties of (i-dimensional spherically symmetric 
solutions. A solution of a d-dimensional Bonnor star is also analyzed in subsection IIII Bl in some 
detail. Its generic properties are shown, as well as its quasi-black hole and its quasi Newtonian 
limits. In Sec. IIVI we present final comments and conclusions. 

II. NEWTON-COULOMB THEORY WITH AN ELECTRICALLY CHARGED FLUID IN 
rf-DIMENSIONAL SPACETIMES (d= n+ 1, n BEING THE NUMBER OF SPACE DIMEN- 
SIONS): WEYL AND MAJUMDAR-PAPAPETROU ANALYSIS, AND BONNOR STAR 
SOLUTIONS 

In (i-dimensional Newtonian gravity coupled to both Coulomb electricity and a charged fluid 
matter, one can find solutions representing charged stars, where here d is the number of spacetime 
dimensions, with d = n + 1, n being the number of space dimensions. The dynamics of such a kind 
of system is governed by the Euler equation, where the gravitational and electric force fields are 
determined conjointly by Newtonian gravity and Coulomb electricity. The fluid can be in static 
equilibrium even with zero pressure and stresses, because the electric repulsion counterbalances 
the gravitational pull if the charge density of the fluid, p e , equals its mass density p m in appro- 
priate units, i.e., \JGd_ p m = ep e , where Gj is Newton's gravitational constant in d dimensions (see 
appendix [A]) , and e = ±1. This condition makes possible to build a distribution of charged dust 
with any shape in neutral equilibrium. Charged fluids with y/Gj p m = ep e are called extremal 
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charge dust fluids. By introducing a convenient boundary one turns the solutions into stars. In 
this section we study some properties of these objects. One can also put some form of pressure, 
either positive or negative, into these systems and find solutions which of course do not obey the 
extremal condition. Solutions with pressure will not be considered here. 

A. The gravitating Newtonian charged dust fluid, and Weyl and Majumdar-Papapetrou 
type analysis 

1. The gravitating Newtonian charged dust fluid 

We consider first the dynamics of a gravitating Newtonian charged fluid in a n = (d — 1)- 
dimensional Euclidean space according to the Euler description. A dust fluid is completely specified 
by its velocity vector, with components Vi, with i = 1, ...,d — 1 (Latin indices run through 1 to 
n = d — 1), and its matter density p m , all being functions of the position vector represented by 
spatial coordinates r,, and of the universal time t. Thus, Vi = Vi(rj,t) and p m = p m (rj,t). The 
basic equations governing the flow of a Newtonian fluid are the continuity equation, which expresses 
mass conservation, and the Euler equation, which expresses momentum conservation. Consider a 
fluid element with mass dm = p m dV, in the {d — l)-dimensional space, dV being the (d — 1)- 
dimensional space volume element. Then, the continuity and the Euler equations may be written 
as 



respectively, where d/dt = d/dt + v % Vi is the convective temporal derivative, Vj is the (d — 1)- 
dimensional gradient operator, Fj is the volumetric external force acting upon the fluid element, 
and the sum convention on indices is adopted. The Newtonian systems we are interested in here 
are gravitating charged fluids distributions in static equilibrium. The fluid is then allowed to have 
some net electric charge, so that the charge of a fluid element is dq = p e dV, p c being the electric 
charge density of the fluid. Thus, there are two independent forces acting on a fluid element, the 
gravitational and electrostatic forces. Both these forces may be derived from scalar potentials, V 
and (f>, respectively, such that one has 



dp, 



+ V, (fay*) =0, 



(1) 



dt 




(2) 



Fi = -puNi V - p c Vi <j) . 



(3) 
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The gravitational potential V is related to the mass density p m by 

V 2 V = S d _ 2 G dPm , (4) 
while the electric potential 4> is related to the charge density p e by 

V 2 <P = S d . 2 Pe, (5) 

where the operator V 2 is the Laplace operator in d — 1 space dimensions, S d -2 is the area of the 
unit sphere in (d — l)-dimensional space given by = 2vr( d - 1 )/ 2 /r(((i - l)/2), r is the usual 

gamma function, and G d is Newton's gravitational constant in d = n + 1 dimensions (see Appendix 
|A]for the definition of Gd)- S d -2 reduces to 4ir in four spacetime dimensions and Eqs. (j3J) and © 
are the natural generalizations of the corresponding three-dimensional Poisson equations for the 
potentials V and 4> to {d — l)-dimensional space. 

We will consider only static systems, so all quantities are functions of the d—1 space coordinates 
only, and the fluid's velocity can be put equal to zero, V{ = 0. Then the Euler equation ((2J) for the 
charged fluid in static equilibrium reads 

p m ViV + p e Vi<j) = 0. (6) 

Equations ar e the important equations for the problem. Eqs. (j3J) and ([S} are the field 

equations that determine the gravitational and the electric potentials once the mass and charge 
densities are given, while Eq. flS} is the equilibrium equation for the system. 



2. Weyl and Majumdar-Papapetrou type analysis 



In vacuum, doing for Newtonian gravity what Weyl did for general relativity 24|, assume now 



an ansatz, i.e., a functional relation, between the gravitational and the electric potential, 

V = V{4>) . (7) 

Eq. ((7J) is the Weyl ansatz which implies that V and eft have the same equipotential surfaces. With 
this ansatz, Weyl originally worked out the Einstein-Maxwell vacuum equations that would follow 
and found that the relativistic potential is a quadratic function of the electric potential. Doing 
the same here in Newtonian gravity, we find that the ansatz ([7]), in vacuum, p m = and p e = 0, 
when put into Eqs. flU) and ©, gives the following equation, (V'f V 2 cf> + V'V" (Vj0) 2 = 0, 
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where the prime stands for the derivative with respect to <f). Thus, since V 2 <fi = in vacuum, and 
(Vj 4>) 2 7^ 0, it follows that , V" = 0, i.e., V{4>) = a o + const x <fi , where ao is an arbitrary constant, 
that without loss of generality can be put to zero. In addition, with our choice of units one has 
that const = —t\[Gd- Thus, 



V(^) = ao-ey / G d ^. (8) 



25] 



This is the Weyl relation for the Newton-Coulomb theory in vacuum. Following Majumdar 
and Papapetrou [3] lead in general relativity, it is interesting to investigate the consequences of 
the linear relation between electric and Newtonian potentials, V = ao + const x (f>, and see what 
happens in presence of matter. Such a relation is a particular case of the Weyl's general ansatz ([7]), 
and is the same as in Eq. (jHJ), i.e., is the same as Weyl's relation for vacuum. It is remarkable that 
the Majumdar-Papapetrou relation is equivalent to Weyl's relation in Newton-Coulomb theory in 
vacuum, while it is not so in general relativity. We use this relation ® to treat also solutions with 

n 

matter, as has been done in general relativity 



In matter, we work out the basic equations using the general form ([8]) and so generalize to 
higher dimensions the analysis in four dimensions done by Bonnor . We show for d — 1 space 
dimensions the interesting result that the equality of mass and charge densities follows from the 
field equations as long as there are no singularities within the charged matter distribution (see 



Bonnor [8|] for Newtonian systems, and De and Raychaudhuri 30(] for general relativistic systems 
in four dimensions). The basic equations ©, and ©, can be rewritten by taking the Weyl 
ansatz (|7|) into account. To begin with, it is convenient to consider first Eq. ([6]), which now reads 
(Pm V + p e ) Vj0 = . So, the two fields V and cfi have the same equipotential surfaces. Since we 
consider Vj <j> 7^ 0, Eq. ([6]) is then equivalent to p m V'+p e = , where again the prime stands for the 
derivative with respect to 4>. By substituting p m from the previous equation into Eq. (J3|), one finds 
(V') 2 V 2 4> + V V" (Vj (ft) 2 = —Sd-2 Pe , where we made use of the assumption V = V{4>). Then, 
with the help of Eq. © one finds Vj [yfz V>) = , where Z is defined as Z = G d - V' 2 . Now, in 
order to have a nonsingular solution with closed boundary it is required that Z = 0, or equivalently, 
iy') 2 = Gd- All equilibrium solutions with (V) 2 7^ Gd with a closed equipotential hypersurface 
S have a singularity within S. The most physically interesting solutions are then those for which 
iV') 2 = Gd- Therefore, the Majumdar-Papapetrou relation for the Newton-Coulomb theory with 
a charged dust fluid is V{4>) = «o — £\[G~d <P ■> where e = ±1, implying, after Eq. ([5D, that the const 
appearing before the potential cj> is indeed —e^JGd- Thus, for the Weyl relation or Majumdar- 
Papapetrou relation (they are the same here), Eq. ©, with the equations p m V + p c = and 
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V(4>) = ao — e^fGd 4> derived above, gives p e = e\JG d p m . This last equation is the Majumdar- 
Papapetrou condition in Newtonian gravity. Observe that the relation between the potentials 
we call Majumdar-Papapetrou relation and the relation between the densities we call Majumdar- 
Papapetrou condition. In the case the Majumdar-Papapetrou condition holds, distributions of 
charged dust of any shape can be put in equilibrium. All the quantities can now be given. Once 
the mass density p m is given, the gravitational potential is determined by the Poisson equation 
([9]) and all the other quantities, including the electro-gravitational Newtonian spacetime structure, 
and possible singularity structure, follow from V and p m . The resulting system of equations can 
be put in the form 

V 2 V = S d ^ 2 G dPm , (9) 



€ 



V, (10) 



vGd 

t\[Gd~p m - (11) 



where the zero points of the potentials are suitably chosen. 



B. Spherical d spacetime (n space) dimensional Newtonian Bonnor star solutions 



1. Equations in spherical coordinates 

We now assume the mass distribution is spherically symmetric, in which case all the dynamical 
variables and fields depend only on the radial coordinate r in (d — l)-dimensional space. Our 
interest here is in spherical solutions to equations ([9])- (jlip . First we define the mass m(r) and the 
electric charge q(r) inside a sphere of radius r, respectively, as 

m(r) = S d _ 2 p m (r)r d - 2 dr, (12) 
Jo 

q(r) = S d - 2 ^ Pc (r)r d - 2 dr. (13) 
Jo 

Equations ([9"j)- (|lip are then conveniently written explicitly in terms of the radial coordinate as 

dV(r) m(r) 

— = Gd ^> (14) 
g(r) or 4>(r) = ^=V(r), (15) 



dr r d ~ 2 ' ^ v ; y/G~ d 

q{r) = ey/Gd~m(r) , (16) 
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where the zero points of the potentials were suitably chosen. In Eq. (|14p there is also a term 
Co/r d ~ 2 which we have put to zero, without loss of generality, i.e., Co = 0. This term can be 
included in the term Gdm{r)/r d ~ 2 by an appropriate choice of the function m(r). 

2. Solutions 

(a) Electrovacuum solutions in d = n + 1 spacetime dimensions 

The solution to Eq. (|14p in vacuum is 

y = 1 G d M 

d-3 r d ~ 3 
M = constant , 

with M representing the Newtonian mass of the source. To complete the solution one must give 
the electric potential <fi, which is obtained from Eq. (|15p . <fi = (d — 2>)^ 1 Q/r d ^ i , with Q being 
the total charge of the source, which in turn satisfies Eq. (fT6|) . Q = e^JUdM. These two equations, 
together with Eqs. (|17p and (|18p . form the set of equations corresponding to the solution of a 
Newtonian Majumdar-Papapetrou vacuum system in n = d — 1 space dimensions. Such a set of 
solutions also follows from the Poisson equations in which the mass and charge densities are Dirac 
delta functions, p m {r) = M 6(r) and p e {r) = QS(r), and Q = e\/UJ M. 

(b) Newtonian Bonnor star solutions in d = n + 1 spacetime dimensions (n space dimensions) 

Now we find a class of solutions to the Newton-Coulomb system with electrically charged fluid 
matter, considering the Majumdar-Papapetrou relation ([8j) (which in the Newtonian case is also 
Weyl's relation) and a spherically symmetric distribution of matter. Upon joining this class of 
solutions to an external vacuum we obtain Bonnor stars in the Newtonian theory. The relevant 
equations are the ones presented in system (fT2|) - (fl6|) . 

Let us call ro, the radius of the star. Physical conditions require the mass density to be a 
continuous function with a finite value at the center of the star. One can choose a mass density 
function p(r) and the remaining functions are then obtained by integrating the appropriate equa- 
tions. For instance, one can give po (r/ro) a , for < r < ro, and make it zero in all the exterior 
region for r > tq. Integration of the Poisson equation (|14p gives a power law function for the 
potential, a\ r a+2 + a,2/r d ~ 2 + 0,3, where the constant 02 is made equal to zero in order to avoid 

21 



(17) 
(18) 



a singularity at r = 0, and the constant 03 is fixed by the matching conditions at r = tq. Alter- 
natively, instead of giving p m (r), one can choose a potential V(r) satisfying reasonable boundary 
conditions, and then obtaining the other functions from it. This is the simplest route, the one 
we follow here. We can choose the following interesting potential, for the interior V{(r), given by 
V\{r) = Co + c\ (r/ro) a + C2 (r/ro)' 3 , for r < ro, where a and /3 are arbitrary constant parameters, 
possibly satisfying some restrictions. The other constants, cq, c\ and C2, are fixed by imposing 
appropriate matching conditions at the surface of the star, r = tq. One can impose that the poten- 
tials are C 1 functions at ro, which is usually done in order to simplify the calculations. This means 
continuity of the potential and continuity of the gravitational field strength. Then, in this case, the 
density has a finite discontinuity at the boundary, falling from some finite value just inside matter 
to zero just outside. Also through the junction conditions above, one can find the constants cq, c\, 
and C2, with one of them arbitrary, c\ say. Here we want to go a step further and impose that the 
potentials are C 2 functions at ro, i.e., continuity of the potential, continuity of its first derivative, 
and continuity of its second derivative. Continuity of the first derivative of the potential means 
that the gravitational field strength is continuous, and continuity of the second derivative means 
that the mass density at the surface of the star is continuous with zero value. Continuity of the 
potential gives Vi(r ) = V c (r ) = -(d - 3)- 1 G d M/r d ~ 3 , where V e (r) = -{d - 2>)~ l G d M/r d ^ z is 
the Newtonian potential in the exterior region, r > ro, with V e zero at infinity, and M is the total 
mass of the star. Continuity of the gravitational field strength gives K'(ru) = V^'(ro) = G d M Jr d ~ 2 . 
Continuity of the mass density at the surface of the star gives j» m (ro) = . With these choices, the 
spherical Newtonian star is described by the following potential 



V 



C M f l+ (d_3)G9 + d-3) 



d - 3 r d - 3 



a(P — a) 



1 



(d-3)(a + d-3) 
P(P ~ a) 
G d M 



1- - 



d - 3 r d ~ 3 ' 
and by the following mass density 

(a + d-3)(/3 + d-3) 



M 



Pn 



(d-3)(/3-a) 



S, 



d-2T0 



d-1 



a-2 



(3-2 



0. 



r < r Q , 
r > r , 

r < r , 
r > ro . 



(19) 



(20) 



In these equations M is the mass and ro is the radius of the star, with M being obtained from Eq. 
(|12p with r = ro, i. e., M = m(ro), and a and (3 are arbitrary constant parameters satisfying the 
restrictions a > 2 and (5 > 2. In addition, the parameters a and (5 must be different from each 



22 



other, a ^ (3, in order that the mass density be finite at r = 0, and the other functions that follow 
from it be also finite there. Such conditions ensure also the positivity of the mass density. Note 
that the quantity M/(Sd-2 r o d ~ 1 ) appears naturally, indeed in Newtonian theory one can define 
the mean density of the matter by p m = (d — 1) Mj (Sd-2 '"o X )- The other quantities, <fi and p e , 
are obtained by substituting the expression for the gravitational potential and for the mass density 
given in Eqs. (|19p - (|20p into Eqs. (|10p and (jlip. respectively. So the class of Bonnor stars is defined 
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FIG. 1: The rescaled Newtonian potential V(r) + 1 as a function of r/fi, where /i = (GdM/(d — 3)) 1 ^ d ~ 3 \ 
for four spacetime dimensions, d = n+ 1, d = A (top- left panel), d — 5 (top-right panel), d = 6 (bottom-left 
panel) and d — 7 (bottom-right panel), and for four different values of the parameter a. The lowest, solid, 
curve is for a = 1, the dot-dashed line is for a = 0.7, the dashed line is for a = 0.4, and the dotted line for 
a = 0.1. 



essentially by equations (|T9|) - (|20|) . through the parameters Gj, M, ro, d, a and /?. In the analysis 
of these Bonnor stars, an important parameter appears, the (i-dimensional generalization of the 
mass to radius ratio of the star, 

G d M 



(21) 



d - 3 r d - 3 ' 

It measures how compact is the star, and is a free parameter in the model. Taking M as a fixed 
parameter, different stars are parameterized by different values of a, which means different values 
of the radius ro- There are no constraints on a for Newtonian stars, it can vary from 0, a highly 
dispersed star, to oo, a point mass, i.e., the limiting configuration here is a Newtonian singularity 
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at r = obeying the Majumdar-Papapetrou condition Q = M. As we will see, in the relativistic 
case a cannot be larger than unity (see also 13. fis|). 
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FIG. 2: The normalized Newtonian mass density p m (T)/X m as a function of r/p, where A m = -rzrPmi Pm 
being the average density (see text), and p = (GdM/(d — 3)) 1 ^ d ^ 3 \ for four different spacetime dimensions 
(d = n + 1 = 4, 5, 6, 7, as indicated) and for four different values of the parameter a.. The solid line is for 
a = 1, the dot-dashed curve is for a = 0.7, the dashed line is for a = 0.4, and the (lowest) dotted line is for 
a = 0.1. The normalized Newtonian mass density p m ('')/A m goes to zero at the surface of the star, defining 
thus the radius tq in each plotted case. 

The relevant functions V(r) + 1, p m (r), </>(r), and p e ( r ), given in terms of the coordinate 
r follow from the above relations. They are dependent on the variable r, and also depend on 
two other arbitrary parameters, the mass and the radius of the star, M and ro, respectively. 
Instead of writing the explicit form of such functions, it is more convenient to plot them for several 
choices of parameters. In the calculations we normalized the coordinate r to the mass parameter 
p, = {Gd M/(d — S)) 1 /^ -3 ) which was kept fixed. In fact, the important parameter to this end 
is the mass to radius ratio a, given by Eq. (I2ip . The function V{r) + 1: The behavior of the 
rescaled potential V(r) + 1 as a function of the rescaled coordinate r/p, for four different values 
of a (a = 0.1, a = 0.4, a = 0.7, and a = 1), and in four different spacetime dimensions (d = 4, 
5, 6, 7) is shown in Fig. [H We plot the rescaled function V(r) + 1 instead of simply V{r) for 
direct comparison with the relativistic case studied later. Now, the parameters a and (3 in the 
solution (fT9j) - ([20j) are free parameters. We have chosen them so that (3 = 3a/2 = 3(d — 3), which 
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is a convenient choice when one studies the counterparts of these solutions in general relativity. 
With this choice, the form of the curves depends on the number of spacetime dimensions d and 
on the parameter a alone. Note that all the interior functions Vj(r, a) match the exterior solutionn 
V c (r) + 1 = 1 — GdM/((d — 3)r d ~ 3 ), each one at a different value of tq. The reason for that is 
because the change of a is made by keeping the mass of the star fixed, while ro varies accordingly. 
The function /0 m (r): Another quantity of interest is the mass density p m (^)- in Fig- El we plot 
p m (r)/A m as a function of the normalized radial coordinate r/p. The density A m is defined as 

Am = [d- 1)7^-3) t/3-a) P m ' wnere tne mean density p m is given by p m = (d - l)M/(S d _ 2 r^' 1 )- For 
our choice of parameters, (3 = 3a/2 = 3(d — 3), one has A m = xrrP m . It is seen that p m (r) is finite 
at r = 0. In fact, with our choice, p m vanishes at r = for all d > 4. In addition it goes to zero 
at the surface of the star, defining thus the radius ro in each plotted case. The behavior of the 
potential <p{r) is simply given by 4>{r) = — (e/y/Gd) V(r), and it is not plotted. The behavior of 
the charge density is p e (r) = £y/G~2p m , and it is not plotted. Note that the potentials, V and 4>, 
are C 2 functions of r, so that the corresponding field strengths are continuous (C 1 functions, in 
fact) through the surface of the star. The mass and charge densities, p m and p e , are C° functions 
vanishing at r = ro. When ro — > one obtains a point charge with a central Newtonian, mild, 
singularity. It is mild because it is not a nasty spacetime singularity, it is a matter singularity only. 

III. EINSTEIN-MAXWELL THEORY WITH AN ELECTRICALLY CHARGED FLUID 
IN d SPACETIME DIMENSIONS (d = n + l): WEYL AND MAJUMDAR-PAPAPETROU 
ANALYSIS, AND BONNOR STAR SOLUTIONS 

In (i-dimensional general relativity coupled to both Maxwell electromagnetism and a charged 
fluid matter one can also find solutions representing charged stars. The fluid can be in static 
relativistic equilibrium if it is made of extremal matter, where the electric repulsion from the 
charge density of the fluid, p e , counterbalances the gravitational pull from its mass density, p m , 
in appropriate units, i.e., y/Gj p m = ep e . Thus, relativistic Bonnor stars in d dimensions can 
also be constructed. Within general relativity the behavior and properties of these solutions is 
much richer, allowing the possibility of quasi-black hole behavior, for a sufficient compact object, 
rather than the point like dull singularity of Newtonian objects. In this section we study some 
properties of relativistic charged fluids in the context of a Majumdar-Papapetrou analysis, and the 
corresponding Bonnor stars. 
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A. The gravitating relativistic charged dust fluid, and Weyl and Majumdar-Papapetrou 
analysis 



1. The relativistic gravitating charged dust fluid 



With the aim of finding exact solutions for d-dimensional Bonnor stars we firstly write the basic 
equations for the Majumdar-Papapetrou systems and analyze their general properties in brief. 
In the following sections we particularize for spherically symmetric spacetimes, show a particular 
solution and study it in some detail. 

The general relativistic analog of the Newtonian charged fluid discussed in the preceding section 
was considered in 



67(. Such a relativistic system is described by the d-dimensional Einstein- 
Maxwell with an electrically charged fluid system of equations which read (we use units such that 
c = l), 

Gfj,u = ^ _ ^ Sd-iGa (E^ + T^ u ) , (22) 
V,F^ = S d . 2 J» , (23) 

with G^u, being the Einstein tensor, such that G^ = — ^g^ u R, R^ being the Ricci tensor and 
R the Ricci scalar. The right hand side of Eq. (|22p bears a universal constant Gd, which in four 
dimensions corresponds to the Newton's gravitational constant (see Appendix [X] for the definition 
of G d ). S d - 2 = 2^- 1 )/ 2 /r(((i- l)/2), where V is the usual gamma function, and the whole factor 
[d — 2) GdSd-2/ {d — 3) corresponds to the 8ttG term in four dimensions. The electromagnetic 
energy-momentum tensor, E^ v , is given by 

where = V ^A v — V U A^, is the electromagnetic gauge field, with V M being the covariant 
derivative. J M , in Eq. ([23]) . is the current density 

J H = PeU^ , (25) 

where p e is the charge density and is the velocity of the fluid in the (i-dimensional spacetime 
with g^u^vy = — 1. Finally, T^ u is the the matter energy- momentum tensor for dust given by 

T^u = PmU^Uv , (26) 

with p m being the energy density of the fluid. In all the above definitions, Greek indexes /i, v, etc., 
run from to d — 1, where represents the time, and the other d — 1 coordinates are spacelike. 
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It is assumed the spacetime is static, in which case the metric can be written in the form 

ds 2 = -W 2 dt 2 + —s-hijC^daP , (27) 
W d ~ :i 

where Latin indices run from 1 to d — 1, hij is the metric in (d — l)-dimensional space, and If is a 
function of the spacelike coordinates x l only. The four-velocity and the gauge field are then given 
respectively by 

u^ = W8l, (28) 

and 

A M = (29) 

where the electric potential ip is also a function of the space coordinates alone. (Note that in the 
definition of we have put a minus sign in front of (p. Although not the usual choice, this is the 
useful choice to compare with the Newtonian case.) 

From the Einstein-Maxwell with a charged dust fluid equations one obtains the following equa- 



tions for W and (p 

V 2 W - (V^) 2 = ^ (V^) 2 + 5 d _ 2 G d W^ Pm , (30) 
VV = 2 ^ V V - S d _ 2 W^ Pe , (31) 
where Vj stands for the covariant derivative with respect to the space metric hij. Making now the 



connection to the Newton-Coulomb theory with a charged dust fluid, one may say that Einstein- 
Maxwell with a charged dust fluid equations, Eqs. ([30]) and ([3T]). correspond to the Poisson 
equations for the gravitational and electric potentials, Eqs. and ©, respectively. Moreover, 
continuity and Euler equations (HD and @ are, in certain sense, analogous to the relativistic 
conservation equations V ^E^ v + V^T^ = , which in turn also follow from the general relativity 
equations. In the present case one has 

p m ViW + Pe Vi(p = 0, (32) 

for the conservation equation, which has the same form as Eq. ([6]). 

2. Weyl and Majumdar-Papapetrou analysis 

In vacuum, the generalization of the Majumdar-Papapetrou system to d-dimensional spacetimes 
was done in [671 ] and, for completeness, we summarize the main properties of such systems here. 
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Following the lines of that work but changing the strategy in order to compare the present analysis 
to the Newton-Coulomb case of previous sections, we assume there is a Weyl-type functional 
relation between the metric potential W and the relativistic electric potential (p 

W = W(if) . (33) 

This is the relativistic Weyl's ansatz. For the sake of comparison to the Newtonian case, let us 
review here the main consequences of the last equation. For the vacuum case, p m = p e = 0. So, 
Eqs. 421 and ^H) can be combined to yield (V^) 2 (W W" + W' 2 - WG d ) = . Since (V^) 2 ^ 0, 
this equation implies in WW" + W' 2 — WGd = 0, which integrates to W 2 = (ao — eVGd (p) 2 + b , 
where ao and 6o are integration constants. This form of the metric potential W is known as the 
Weyl potential or, in our context, the Weyl relation. Moreover, in the particular case where bo = 0, 
W 2 assumes the form of a perfect square so that 



W = oq - ey/G d <p, (34) 

where e = ±1, and without loss of generality we kept the plus sign when taking the square root 
of W 2 . In general relativity, this form of W is known as the Majumdar-Papapetrou potential, and 
one usually refers to Eq. (|34p as the Majumdar-Papapetrou relation. 



In matter, we now render into d dimensions De and Raychaudhuri's theorem [30j] (see [72| for 
the generalization of it for systems with pressure). To begin with, one substitutes of Eq. (|33[) 
into the conservation equation (|32p and finds (p m W + p c ) Vnp = which, for Vi<p ^ 0, is then 
equivalent to 

p m W + p e = 0, (35) 

where the prime denotes derivative with respect to <p. This is the general relativistic analog to 
the equilibrium equation of Newtonian theory, cf. the equation p m V + p e = derived in Sec. 
Ill A 21 Using Eq. ([35]) . it is also readily seen that, together with Eqs. ([30]) and (f3T|) . it implies 
in Vj (\/ r ZX7 l ip/W) = where here Z = Gd — W' 2 . This equation is to be compared to its 



Newtonian analog and, in fact, has the same form. So it is possible to generalize the theorem by 



De and Raychaudhuri [30j to higher dimensions (see also 



14l|). According to such a theorem, in 



order to have charged dust solutions satisfying Weyl hypothesis without singularities, the quantity 
Z must vanish. This implies W' 2 = Gd, as in the Newton-Coulomb with electric matter theory, so 
that the result is the Majumdar-Papapetrou relation, the same as in the relativistic vacuum case 
(see Eq. (fM)) ). W = ao — ey/G~d~(p, with e = ±1, and ao being an integration constant. Then, 
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substituting W from the latter equation into Eq. (|35|) gives p e = e y/Gj Pm as in the Newtonian 
case. To sum up, let us write the resulting equations for the important functions W, tp, p m , and 
p c . In order to get a field equation similar to Poisson equation ©, it is convenient to introduce a 
new potential U such that 



u = w- (36) 



The relevant equations are then 



O d-l 
7^ 



V 2 U = -S d ^ 2 G d U** Pmi (37) 
V = ^fl-^> (38) 



p e = e \/~Gd~p xa , (39) 

where an arbitrary constant in the potentials was adjusted to unity. Some special solutions to 
these type of systems are going to be analyzed in the next sections. Eq. (|39p is the Majumdar- 
Papapetrou condition. Note, that these equations can be compared to the Newton-Coulomb with 
an electrically charged fluid case. In fact, taking the Newtonian limit in which U — 1 — V, with 
\V\ « 1, one sees that Eqs. §7}, ([38]) and ([39]) reduce exactly to Eqs. ©, dTUJ) and ([TT]) . 
respectively. 



B. Spherical d spacetime dimensional relativistic Bonnor star solutions 

1. Equations in spherical coordinates 

In what follows we confine attention to spherically symmetric static spacetimes and write the 
foregoing equations in isotropic and Schwarzschild spherical coordinates. 

Equations in isotropic coordinates: The starting point is the metric (j27j) . which with Eq. (|36|) now 
reads 

ds 2 = -U~ 2 dt 2 + L/3=5 (dR 2 + R 2 dn 2 d _ 2 ) , (40) 

with d£l 2 d 2 being the metric on the unit [d — 2)-dimensional sphere S d ~ 2 . U is now a function of 
the radial coordinate R only, and it obeys 

" 8 " G "-*" D *- (41> 
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which is obtained from Eq. (|37[) . The matter and charge densities are also functions of R only, 
Pm = Pm(R) and p e = Pc{R), and they are related to each other through Eq. (|3"9"j) . From Eq. ([Hj) . 
and in analogy with the Newtonian theory, define the mass function m{R) and the charge function 
q(R), (see other mass function definitions in Appendix [B]) , as 

m{R) = S d - 2 / p m {R)U(R)^R d ~ 2 dR, (42) 



o 



pR d _ 1 

q(R) = S d _ 2 / Pe (R)U(R)^R d - 2 dR. (43) 
■/ o 



Eqs. (f37]) - (f39j) may then be written as 



dU(R) m(R) 

<K ' (R) - -U(R)-^, or m = -^(l-T^), (45) 



q(R) = e^m{R), (46) 

where arbitrary constants in the potentials were set to one. In addition, terms of the form 
const/it" 1 ' -2 in Eqs. (|44p and (|45p were not written explicitly since they are implicitly absorbed in 
those equations, and moreover they should be put to zero as the fields shall be regular functions 
of the radial coordinate R. Eqs. ()44|) -(ll6l) can then be compared to the Newton-Coulomb with 
an electrically charged fluid case. In fact, taking the Newtonian limit in which U — 1 — V, with 
\V\ « 1, and R ~ r, one sees that Eqs. (jUty - flM]) reduce exactly to Eqs. (fT4"|) - (fl~6]) . respectively. 
For future reference we write here the Kretschmann (/C) and Ricci (1Z) scalars for the metric (|4Up : 

d-l AU" 2 f (3rf-8)(4rf- 11) \ 2U' A _ (rf-2)(2rf-5) AU' 2 U" 

d^~^ + { + (d=W (^-3) 2 

8S d _ 2 Pm ( U' 2 5 d _ 2 d-2 ft\ 



d-3 V U 2 u d-3' 

25 d _ 2 d-4 U' 2 

U = 7^ GdPm -^3^T' (48) 

where the prime stands for the derivative with respect to R. From this it is seen that spacetime 
singularities occur at points where U = 0, as long as the derivatives of U do not vanish at the 
same points as U does. Although the field equations are easily written and solved by working in 
harmonic coordinates, the physical interpretation of the solutions is clearer if one uses Schwarzschild 
coordinates. 

Equations in Schwarzschild coordinates: In Schwarzschild coordinates the line element reads 

ds 2 = -B 2 dt 2 + A 2 dr 2 + r 2 dQ 2 d _ 2 , (49) 
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where B = B{r) and A = A(r), r being the Schwarzschildean radial coordinate. By comparing 
Eq. PU|) to Eq. (|39J) , we see that the radial coordinates in the two systems are related by 

r d - 3 = UR d - 3 , (50) 

and that the metric potentials are related by 

B = ±, (51) 

and 

A-l--?— ^ (52) 

A ~ l d-3Udr- [b2 > 

Eq. (|50p gives r as a function of R. Although this implicitly determines R as a function of r, it is 
only in special cases that this relation can be worked out explicitly. For the sake of completeness, 
we present here the Schwarzschild coordinate form of the field equations. With the metric in the 
form of Eq. ([15]) . Eq. (}3T|) turns into 

\Ur^4- d 4\=S d - 2 G d r^ Pm . (53) 



A dr V AB dr 

This is, in fact, the equation for the potential B, since A is not independent of B. Namely, Eqs. 
fl5Tj) and d52j) give 

a 1 r dB . 

which is a consequence of the Majumdar-Papapetrou condition in a fluid with vanishing stresses. 
At last, the electric functions are expressed in Schwarzschild coordinates. No effort is needed 
to obtain the electric charge density since it is proportional the the mass density. The electric 
potential (p(r) comes after Eqs. ([38]) and ([5T|) . i.e., 

<P = (1 - B) , (55) 

where, as usual, the arbitrary constant was set to unity. Now, defining M{r) and Q(r), (see for 
comparison other mass definitions in Appendix |B]) , as 

M(r) = 5 d _ 2 ^ Pm (r)A(r)r d - 2 dr, (56) 
Jo 

Q(r) = S d _ 2 Pc (r)A(r)r d - 2 dr, (57) 
J o 
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Eqs. (|37j) - (|39|) may be written as 



dB(r) M(r) 

= G d A{r) B(r)^l , (58) 

<lr ' ir] -A(r)B(r)^±, or <p(r) = (1 - B(r)) , (59) 



dr r d 2 yG~d 

Q(r) = ey^M(r), (60) 

and A{r) is given in terms of B(r) by Eq. (|54p . These are the fundamental equations in 
Schwarzschild coordinates. The Newtonian limit is obtained by noticing that for weak gravity 
fields one has that the metric functions B(r) and A(r) are close to unity, B(r) = 1 + 5B(r), and 
A(r) = 1 + 5A(r), with 5 indicating small quantities. Hence, to the first order approximation, the 
above equations reduce respectively to Eqs. (fH| - (fTI)]) . 



2. Solutions 



661 ] . The general 



(a) Electrovacuum solutions in d spacetime dimensions 

As a first example and to set up notation let us report here on the case of d-dimensional 
vacuum Majumdar-Papapetrou solutions. These are nothing but the extreme Reissner-Nordstrom 
spacetimes generalized to higher dimensions that were first studied in Ref. 
solution of Eq. (|44p in vacuum is usually written in the form 

U=1 + —3B^> (61) 
M = constant, (62) 

where M is an integration constant equal to the total mass of the source. The electric potential 
follows from Eq. (|45p . cp = e (1 — 1/U) /y/Gd, and the electric charge is related to the total mass of 
the source by Q = e \J~G~2 M, as required by the Majumdar-Papapetrou condition and in agreement 
with Eq. (|46p . The corresponding spacetime metric is 

Using Eqs. ()50)) - (j52[ ) we find the relation between r and R, given by 

r d-3 = R d-3 + . (64) 

d 3 
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One also finds that B = -4 = [ ^—r | = ( 1 — ^-1^) which leads to the metric for an 

extreme Reissner-Nordstrom black hole with mass and charge equal to M, and holds for all d > 4, 

G d M \ 2 2 , <ir 2 



The coordinate r can be extended up to r = 0, which is in fact a singularity. This is seen from the 
Ricci and Kretschmann scalars which are, respectively, 

; , (d-l)(d-2) 2 G^M 2 / (3d-8)(4d-ll) \ G\M± 

L " ' — — + + ^3 )^) 

(66) 



d - 3 r 2 ( rf -!) V (d-3) 3 J r 4 ( d - 2 ) 

(2d-5)((i-2) 2 G 3 d M 3 
(d - 3) 2 r( M - 5 ) ' 



d - 4 G 2 M 2 

* = -dTs;*w' (67) 

where we used Eqs. (|47p . (|48p . (|50p and (|64p . The region of the spacetime which in Schwarzschild 
coordinates corresponds to < r < (GdM/(d — S)) 1 /^ -3 ) is not covered by the isotropic coordi- 
nates. The maximal analytical extension of these vacuum solutions representing extremal black 
holes can then be found following the usual methods. 

(b) Relativistic Bonnor star solutions in d spacetime dimensions 

Interesting exact solutions in the context of Majumdar-Papapetrou relativistic systems are 
the Bonnor stars, see now specifically [(], Q, [3], which are spherically symmetric distributions of 
a charged dust fluid satisfying the Einstein-Maxwell with matter equations in four-dimensional 
spacetimes. The d-dimensional version of such kind of stars are solution to Eq. (|42l) . or Eq. (I53p . 
with appropriate boundary and matching conditions. We look for solutions using the equations 
in harmonic coordinates, and then do the analysis in Schwarzschild coordinates. In order to find 
solutions to Eq. (l4"2j) . a first, possible, procedure is to provide the mass density as a function of 
the radial coordinate, p m = p m (R). This is the procedure usually adopted, because it furnishes by 
construction physically acceptable mass distribution for the star. In the present case, however, such 
a strategy is not advisable because it results in a second order non-linear differential equation for 
U(R), whose solutions can be found just after fixing the number of dimensions of the spacetime. 
A second procedure, of no interest in the Newtonian case, but valuable here, is to choose the 
energy density profile in such a way to transform Eq. (|42|) into an equation whose solutions are 



known, such as the case of sine-Gordon equation used in Ref. 



451 ] . or transforming it into a 



linear equation, so that one can use the well known methods to solve ordinary linear second order 
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differential equations to find solutions. A third alternative procedure is to fix a priori the metric 
potential U = U(R), and then determining the other physical quantities that follow from it. This 
is the strategy we follow here, it allows to write the solutions in closed form, and it is the same 
strategy as the one opted for in the Newtonian Bonnor stars studied above. 



(i) Solutions with smooth boundary conditions and some special solutions 



First we make the analysis in isotropic coordinates. We consider the general relativistic analog 
of the one studied in Sec. Ill B 2f b) (see also [(], tIQ])- We then choose 



U 



Ui = co + ciR a + c 2 bP , R<Ro, 
1 G d M 



(68) 



U, 



1 + 



R > Rq . 



d-3 R d - 3 ' 

where a and (3 are real numbers and -Ro shall be identified as the surface of the star. The arbitrary 
constants cq, c\, and C2 are fixed in such a way to guarantee the matching conditions at the surface 
of the star, R = Rq. Bonnor [(], tIQ] imposed U to be a C 1 function and the energy density to be 
a step function at the boundary. In this case one can verify that the constants are given by cq = 



1 + 



G d M ( P+d-3 



+ ci Rq, ci one can take as arbitrary, and C2 



G d M 



pd R { 







To reproduce Bonnor's choice for U [121 ] one puts d = 4, ci = and j3 = n (where n was the 



letter chosen for the exponent in 



12j|). Of course, if one wishes, one can choose U to be of any 



degree of differentiability at the boundary. Since it is interesting to test whether this choice of 
differentiability has any important influence on the properties of the star one can, still in the spirit 
of Bonnor, go a step further and instead of choosing U as a C 1 function, impose U to be a C 2 
function of R. As a bonus, one gets in addition, that the energy density is a C° function, i.e., 
continuous at the boundary Rq, indeed zero, which is more in accord with the usual properties of 
stars. For a C 2 choice for U there are no free constants and one finds, 



C-2 



1 + 



1 G d M 
d - 3 Rt 3 



1 + 



d-3 /p+d-3 a+d-3 



j3 — a 



a 







{0 + d-3) G d M 
a((3-a) R«+ d - 3 
(a + d-3) G d M 



(69) 
(70) 
(71) 
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It then follows the potentials U{ and U e are 



U 



Ui =1 + 



G n 



M 



d-3Rd~3 



1 + 



(d-3)(/3 + d-3) 
a(/3 — q) 



R 

Rq 



(d-3)(a + d- 3) 



U e =l + 



PW ~ a) 
1 G d M 



R 
Rn 



d-3 R d s ' 
Eq. (|37H then gives the mass density 

'(a + d- 3)(/3 + d-3) M 



(d-3)(/3-a) 



a-2 



R_ 

Rn 



(3-2 



d-1 » 
If d-3 



R — Rn 
R > Rq 

R< R 



(72) 



o 



(73) 



0. 



R > Rq . 



In the region outside the mass distribution, the solution takes the extreme Reissner-Nordstrom 
form (|63p . as expected. Since U is a C 2 function, the spacetime metric satisfies the Israel matching 
conditions at R = i?o- In order that p m be a well defined function and everywhere non- negative we 
must have a,(3>2. The quantity M / (Sd-2 Ro d l ) appears naturally with units of mass density. 
Note that the electric potential <p and the the electric density p e can be found directly from Eqs. 
(I45p and (|46H . The function ip is a continuous C 2 function through the surface of the star, which 
means the field strength is C 1 and the charge density is C°. Moreover, using Eqs. (|42l) and (|73p 
one finds that indeed M = m(Ro), making the whole procedure a consistent one. This Bonnor star 
solution looks like the Newtonian star studied in Sec. Ill B 21 In fact, the resulting mass density, 
Eq. (f73l) . resembles the function given by Eqs. (f20l) . 



Second, we make the analysis in Schwarzschild coordinates. Schwarzschild coordinates are 
interesting to analyze the physical properties of the spherical solutions found above. Eqs. (|50p and 
(|72p establish the relation between the harmonic radial coordinate R and the Schwarzschild radial 
coordinate r 



„d-3 



coR d - 3 + Cl R a+d - 3 + c 2 R p+d ~ 3 , R<R Q , 



R d - 3 + 



1 



GdM , 



R > R 



(74) 



o • 



d-3 

These relations furnish R as a function of r, R = f(r), which is in fact defined by two functions. 
Let us call them respectively /i(r), for the internal region, and f e (f), for the external region. The 
surface of the star, defined by R = Ro, is obtained in terms of the Schwarzschild coordinates, by 
imposing the continuity of the function r(R) through such a surface, i.e., 

1 



d-3 



p d— 3 



Ui{Rn 



t> d—3 
Hq 



U C (R 



r> d—3 
tin 



d-3 



GdM . 



(75) 
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where U{ and U e are denned by Eq. (|72p . The aim now is to find the metric potentials B and A as 
functions of r. In order to do that one needs to find the functions /i(r) and f e (r), which is done 
by solving Eqs. (|74p for R. For r < ro one has 

1 / .„ f p\-i 



B i (r) = T —r=[c + c 1 f i a + c 2 f i p ) , (76) 
Mr) = 1 + t ^ (W" + 0<*/f) f c o + ci/r + c 2 /f) _1 , (77) 



fi dr 

with fi = fi(r) being a suitable solution of the following algebraic equation 

c 2 fi P+d ~ 3 + cif« +d - 3 + c ft 3 ~ r d ~ ?J = . (78) 

The constants Co, c\ and c 2 are now to be written in terms of ro, instead of in terms of Rq. The 
corresponding expressions are obtained by substituting R =(r d_3 -G (i M/(d-3)) 1 /(^ 3 ) into Eqs. 
([69]) - ([7T]) . For r > ro one has 

1 1 G d M , N 

Third, we find some special solutions with a simple algebraic structure. Generally, the only 
way of finding the solutions to Eq. (|78p is by specifying the values of the parameters a and (5 and 
the number of spacetime dimensions d. Even in that case, in general, only numerical solutions are 
possible to find and we do not perform such an analysis here. There are, however, some special 
values of a and for which Eq. (|78p can be solved exactly for fi(r). Thus, in order to investigate 
some more properties of o?-dimensional Bonnor stars, we consider a particular case that can be 
dealt with algebraically. For instance, one may choose 

f3 = - a = 3(d - 3) , (80) 

so that one finds a fourth degree polynomial equation to solve for R d ~ 3 : 

c 2 (R d - 3 ) 4 + ci (R d - 3 ) 3 + c R d - 3 - r d ~ 3 = , (81) 



where now the coefficients cq, c\ and C3 are simplified to 

2G d M 2G d M G d M 

This polynomial equation can be solved in terms of radicals, and the physical quantities can then 
be expressed explicitly in terms of the coordinate r. In order to condense expressions, we first 
define the parameter a by 

G d M 



d — 3 ro 
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( 83 ) 



with < a < 1. As in the case of Newtonian stars (see Eg. (|21|) ) this parameter measures how 
compact is the star and it is useful to parameterize the numerical solutions. Further, we define 



1 / 1\ 2 1 / r \ d ~ 3 , , 1 1 1 / r N 



b(r) = — 1 + - - — - , c(r) = - + ., , 

16 V ay 4a V r o/ 6 6a 3a \ro 



eM = (6(r) + ^(r)] 2 -[c(r)] 3 ) , s(r) = ^1 + 2 e(r) + , 
where we have used the relation Ro d ~ 3 = ro d ~ 3 — t4M. Then, the solution for i?(r) is 



(84) 



d-3 

( 



R(r) d - 3 



r d-s_G<L Mj r>rQ 
a — 3 



(85) 



Fourth, the relevant functions B(r), A(r), Pm(r), <p(r), and p e (r), given in terms of the 
Schwarzschild coordinates follow from the above relations. They are dependent on the variable 
r, and also depend on two other arbitrary parameters, the mass and the radius of the star, M and 
ro, respectively. Instead of writing the explicit form of such functions, which are cumbersome, it 
is more convenient to plot them for several choices of parameters. In the calculations we normal- 
ized the coordinate r to the mass parameter \i = (Gd M/ (d — S)) 1 /^ -3 ) which was kept fixed. In 
fact, the important parameter to this end is the mass to radius ratio a, given by Eq. (|83p . which 
measures how relativistic is the system. Here we have the constraint < a < 1, and for small a 
the system is Newtonian, while for a close to unity it is fully relativistic. The function B{r): The 
simplest function to be found in Schwarzschild coordinates is the metric potential B(r), which is 
immediately obtained through the relation B(r) = l/U(r). Fig. Oshows B(r) as function of r/(i in 
d = 4, 5, 6, 7, as indicated. It is also seen in that figure the behavior of B(r, a) as a function of a, 
for different values of the parameter a, as shown by the four curves in each graph. All the interior 
functions B\(r,a) match the exterior extreme Reissner-Nordstrom solution B e (r) = 1 — (/x/r) a!_3 , 
each one at a different value of rg. The reason for that is because the change of a is made by 
keeping the mass of the star fixed, while ro varies accordingly. Notice also that for a — * 1 the 
function B[(r) approaches zero in the whole region interior to r = ro, meaning that the redshift 
with respect to infinity is infinite. For the extreme value (a = 1) the mass and the charge of the 
charged star are concentrated inside a quasihorizon at r = ro. In this limit, the spacetime solution 
is a quasi-black hole, similar to what was found for four-dimensional spacetimes (see [3]-{l3]). 
There are no singularities inside ro, the curvature is finite, so are the mass and charge densities of 
the charged dust (see also item (ii) below) . It can also be seen the Newtonian limit of the solution 
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B(r) (d=4) B(r) (d=5 ) 




2 4 6 8 0.5 1 1.5 2 2.5 3 fi 



FIG. 3: The metric potential B(r) as a function of r/fi, where /x = (G d M/(d - 3)) 1/(d ~ 3) , for d = 4, 5, 6, 7, 
and for four values of a in each graph (from top to bottom: a = 0.1, a = 0.4, a = 0.7 and a = 1). 



1 1 

Mr) (d=4) A(r) (d=5) 




0.5 1 1.5 2 2.5 3 3.5 4 ^ 0.5 1 1.5 2 2.5 3 ^ 



FIG. 4: The metric potential l/A(r) as a function of r/fi, for d = 4, 5, 6, 7, and for four values of a in each 
graph (from top to bottom: a = 0.1, a = 0.4, a = 0.7 and a = 1). 
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FIG. 5: The normalized relativistic mass density p m (r)/X m as a function of v j '^x, where A m — jiPmi /^m 
being a kind of average density (see text), and /i = (GdM/(d — 3)Y^ d ~ 3 \ for the cases d = 4, 5, 6, 7 (as 
indicated), and with a = 1 (upper curve), a = 0.7 (dot-dashed curve), a = 0.4 (dashed curve), and a = 0.1 
(lowest curve) for each d. The normalized relativistic mass density /9 m (r)/A m goes to zero at the surface of 
the star, defining thus the radius ro in each plotted case. 



by comparing the curves for the smaller values of a in Fig. [3] with the corresponding curves for the 
Newtonian potential, Fig. [1] (see also item (iii) below). The function A{r): The behavior of the 
other metric potential A(r) is seen in Fig. HI where we plot 1/A against r / ' [i for the same values 
of d and a as in Fig. [3j The quasi-black hole formation is seen in this case as 1/A{r) going to zero 
at r = ro for a — > 1. It appears in the figure as the sharp elbow in the solid line (lowest) curve 
showed in the graph. The exterior function is A e {r) = 1/B c (r), and all the inner functions A[(r, a) 
for different a match A c {r) at a particular value of ro. The function p m (r): Another quantity of 
interest is the mass density p m (r)- In Fig. [5] we plot the normalized mass density p m {r)/\ m as a 
function of the normalized radial coordinate r/fi. Here A m is defined as A m = rarw^r^mzf] Pm, 
where p m , a kind of average density, here is given by p m = (d — 1)M / '{Sd-2 r o _1 )- F° r our choice 
of parameters, see Eq. (f80j) . one has A m = -j^j p m . We plot p m (r)/X m against r/p for the same 
values of d and a as in Figs. [3] and HI Notice that for d > 4 the general properties of this func- 
tion do not depend upon the specific value of d. It is clearly seen that p m ( r ) is finite at r = 0. 
In fact, with our choice, p m vanishes at r = for all d > 4. In addition it goes to zero at the 
surface of the star, defining thus the radius ro in each plotted case. Moreover, the mass density 
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is everywhere well defined even in the quasi-black hole limit. The comparison to the Newtonian 
case can be done considering the curves for small a in Fig. and comparing the corresponding 
curves in Fig. [2j (see below item (hi)). The functions ip(r) and p e { r )' The other two functions, the 
electric potential tp(r) and the electric charge density p e (r), are so close related to the respective 
gravitational quantities B{r) and p m (r), that no plot need to be drawn for them. In fact, they are 
promptly obtained from their relations to the functions studied above (see Eqs. (|55p and (|39p ). 
namely, (p(r) = e(B(r) - l)/y/UJ, and p e (r) = e^/G d ~p m (r) . 

(ii) The quasi-black hole limit 

For the full relativistic limit, a = 1 — e, with e « 1, it is clear from the previous plots that 
the function 1/A(r) attains a minimum at r/p = 1 + e, such that 1/A(r) = e, where again, 
p = (Gd M/(d — S)) 1 /^ -3 ) . Also, for such a small but non-zero e the configuration is regular 
everywhere with a non- vanishing metric function B. Moreover, in the limit e — > the interior 
metric potential B{ obeys, B[ — > 0, for all r/p < 1. These three features define a quasi-black 
hole, see [it] . \v\ . These three features imply, among other things, that (a) there are infinite 
redshift whole regions rather than surfaces, (b) the object displays naked behavior, i.e., generation 
of infinite tidal forces in a freely falling frame, (c) outer and inner regions become impenetrable 
and disjoint, and (d) for external distant observers the spacetime is indistinguishable from that 
of extremal black holes. The quasi-black hole is on the verge of forming an event horizon, but it 
never forms one, instead, a quasihorizon appears. 

It is of interest to see that in the quasi-black hole limit the metric is well defined and everywhere 
regular. We check this for the interior. One defines, from the isotropic radial coordinate R, a new 
spatial coordinate x by 

R 

x = — , < x < 1 , (86) 

Ro 

from which one sees that the surface of the star is now located at x = 1. Substituting this 
transformation into the interior metric functions and choosing a new time coordinate T according 
to 

the interior metric is now 

ds 2 = _U-2 dT 2 + ( GdMY^ 3) uV(d-3) ( dx 2 + ^2^) (8g) 



d 
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where 



U = 1 + (d - 3) 



(89) 



a(/3-a) ^ " y ' P(a-0) 

This metric is regular throughout the interior region and also at the surface of the star. Moreover, 
at x = 1 one has U = 1. This means that even in the quasi-black hole limit the surface of the 
star is timelike for internal observers. On the other hand, one can verify that, being the exterior 
metric the extremal Reissner-Nordstrom metric, the quasi-black hole limit gives a null surface for 
external observers. There is thus a mismatch, implying in this case that the interior and exterior 



(Hi) The quasi- Newtonian limit: the Newtonian Bonnor stars discussed previously 

It is expected that in the weak field approximation a relativistic Bonnor star reduces to a 
Newtonian Bonnor star. Here, we show that indeed the relativistic star studied in this section, i.e., 
Sec. HUB 2\ reduces to the Newtonian star studied in Sec. Ill B 21 

In the relativistic theory two coordinate systems are involved in the solutions, the isotropic and 
the Schwarzschild spherical coordinates. Initially we show that to first order approximation in the 
weak field limit the two coordinate systems are identical. In order to deal with the issue we take 
the special case considered in paragraph (b)(i) of subsection IIII B 2l The weak field limit inside 
the spherical star corresponds to small values of the parameter a = G^M/((d — 3)ro d ~ 3 ). Hence, 
considering the approximation of Eq. (|85|) up to the first order in o it follows 



R 



d-3 



jyd-3 



r d—3 
d~3 



1 



G d M 



d-3 r^~ 3 



r d— 3 
d-3 



+ 



„d— 3 

d-3 



(90) 



At the lowest order approximation it results in 

R d-3 



d-3 



r d—3 
d-3 



(91) 



as expected. Therefore, when comparing the first order approximation of the relativistic solution 
to the Newtonian solution one may work with the isotropic coordinates, identifying the radial 
coordinate R with the Newtonian radial coordinate r. 

The next step is obtaining the potentials and the densities in the weak field approximation and 
comparing them to the Newtonian case. In such a limit one has the relation U = 1 — V, where V 
is the Newtonian potential. Now using the relation (I91j) and Eq. (I72p one can write Ui up to first 
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order in M/ro 



d-3 



1 + _3lJL (i + (d-3)(/3 + d- 3) 



d-3rt 3 



a{(3 — a) 
(d-3)(a + d-3) 



1 



r 
ro 
r 



+ 



(92) 



P(J3 - a) 

From this equation, and from the exterior solution U e , one then finds the potential in Newtonian 



approximation, V = 1 — U, as 



V 



Vi 



M 



d-3 r^" 3 



(d - 3)((3 + d - 3) 



a ((3 — a) 



1-1^ 



+ 



(d-3)(a + d-3) 
/?(/?- a) 
G d M 



1-1^ 



r > ro 



(93) 



d-3 r rf " 3 ' 

The resulting expression is to be compared to the gravitational potential of the Newtonian star 
as given in Eq. (|19p . The two expression become identical if one identifies the gravitational 
constant Gd, the radial coordinate r, and the mass of the star M in both equations. We have 
already shown that, in the weak field approximation, it results R = r + 0(M /r -3 ), and also 
U(R) = U(r) = l + 0(M/rQ~ 3 ). Therefore, substituting such results into Eq. (fT3"j) we find the first 



order approximation for the relativistic mass density. 

' (a + d-3){(3 + d-3) M 



(d-3)(J3-a) 



Of-2 



P-2 



0. 



r < ro , 



r > ro 



(94) 



In order for this result to be identical to Eq. (|20p the mass M and the coordinate r must be 
the same in both equations. It is then straightforward to show that the weak field limits of other 
relativistic quantities such as the metric functions B{r) and A(r), the electric charge density, and 
electric potential all agree with their Newtonian counterparts, as expected. 

Notice that units have been normalized in such a way that the gravitational coupling constant 
in Einstein equations equals to the Newtonian gravitational coupling constant in Poisson equation 
(see Appendix [A]) . Furthermore, the mass densities carry identical units and normalizations due 
to the similarity between Poisson equation for Newtonian gravity, Eq. (JH), and the corresponding 
equation coming from Majumdar-Papapetrou relativistic system, Eq. (|37|). 
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IV. CONCLUSIONS 



We have studied d-dimensional Bonnor star solutions, spherical distributions of extremal charged 
dust joined to extremal charged vacua, both in Newtonian gravity and general relativity. We have 
found that the relativistic solutions present many interesting properties such as forming an extreme 
d-dimensional quasi-black hole, when the mass to radius ratio reaches a critical value. We have also 
found that the Newtonian solutions are limiting cases of the relativistic ones. In this connection it 
is interesting to note that the Bonnor star solutions in Majumdar-Papapetrou Newtonian gravity, 
when contrasted to those Bonnor solutions in Majumdar-Papapetrou general relativity, display 
clearly the departing of the high density structures that may arise in the strong field regime of 
each theory, mild singularities in one theory, quasi-black holes in the other. Moreover, whereas 
there are no solutions for Newtonian stars supported by degenerate pressure in higher dimensions, 
and so no general relativistic solutions also, higher dimensional Bonnor stars, supported by electric 
repulsion do indeed have solutions. This means that the existence of stars in higher dimensions 
depends on the number of dimensions itself, and on the underlying field content of those stars, as 
expected. 
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APPENDIX A: NEWTON'S GRAVITATIONAL CONSTANT G d IN d SPACETIME DI- 
MENSIONS 



Within Newtonian gravity, the Poisson equation for the gravitational field is given by 

V 2 V = k Pm , (Al) 

where k is a constant, related to Newton's gravitational constant Gd in d spacetime dimensions, 
to be determined. Integrating over the space volume V and using the Gauss theorem, one obtains 

I V 2 Vd d ~ 1 x= I V i Vn i dS d _ 2 = k [ p m d d ~ 1 x = kM , (A2) 

where Sd-2 is the boundary surface surrounding the volume V, and n l is the unit normal to the 
surface Sd-2- Considering now spherical symmetry, i.e., 

ViVn 1 = -g r , (A3) 

where g r is defined to be the radial component of the gravitational field, one finds 

ViV rjdS d ^ = -g r S d -2 r d ' 2 . (A4) 



(A5) 



' Sd-2 

Then (TX2]) and flU]) yield 

k M 



Sd-2 r d ~ 2 ' 



The choice in 



59] for k is given by 

k = GdSd-2 ■ (A6) 
This is an interesting choice because it gives 

G d M 

9r = ~^d~2 > ( A7 ) 

i.e., a straight generalization of Newton's force law to d spacetime dimensions, although it puts 
Einstein's equation into a slightly awkward form, 

G^u = ^ _ g Sd-2 Gd T^v , (A8) 



68] for 



where is the Einstein tensor and T^ u is the energy-momentum tensor. The choice in 
k is given by 

k = ^G d d i . (A9) 
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This is also interesting choice because although it gives 

8vG d d-3 M 

9r = -^—2—*> (A10) 

Einstein's equation are written as 

Gp,v = 8vr G d T^ v , (All) 

i.e., a straight generalization of Einstein's equation to d spacetime dimensions. Both definitions of 
k give the correct definition in four dimensions for G d= 4 = G4 = G. In this paper we have opted 
for the definition ([A"6|) . which yields (TAT|) and §M\\ . 



APPENDIX B: MASS DEFINITIONS 

1. Mass functions in isotropic coordinates 

Throughout the paper we used the mass function m{R) defined in Eq. (|42p . In the literature it 
is sometimes defined another mass function M(R). The connection between the two definitions is 
given below. Using Eq. (|4"2j) one gets 

17(12) = l- Gd j R ^ dR, (Bl) 

where an integration constant has been made equal to unity. Eq. (|Bip is consistent with the usual 
form of the potential U outside the mass and charge distributions, i.e., R > Rq. In fact, if we take 
m(R) = M = constant, Eq. ([HI]) yields U(R) = 1 + G d M/((d - 3)R d ~ 3 ), where M is total mass 
of the source. The other mass function of a charged dust distribution M(R) can then be defined 
in analogy with the result for vacuum. This is done by taking U(R) inside the dust in the same 
form as outside, 

, . G d M(R) , . 

Hence, it follows the relation 

M(R) = -(d-3)R d - 3 J R ^^-dR. (B3) 

And so, one sees that the two masses m{R) and M(R) are in general different from each other. The 
two definitions agree just in the region outside the dust fluid, in which case m(R) = M{R) = M is 
the total mass of the source. 
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2. Mass functions in Schwarzschild coordinates 



The mass definition in Schwarzschild coordinates used in the paper is given by Eq. (|56p . Besides 
such a definition, there is a different route to define another mass function M(r). Usually, in the 
literature the mass within a certain sphere of radius r, M(r), is defined through a relation of the 
form 

' (B4) 



1 _ G d MM 
1 d-3~FX=Z 



Interestingly, this mass coincides with Ai(r) as defined in Eq. (|56[) . This can be shown as follows. 
From the last equation it follows 

M{r) = d —^r d -*(l-\\ . (B5) 



G d V A, 

Moreover, using the expression for dB/dr in terms of A obtained from (|54h one gets 

Therefore, comparing Eqs. (|B5|) and (|B6P one obtains 

r d ~ 2 dB 

^ = ^M(r). (B7) 

Substituting this result into Eq. (|53|) and integrating one has 

M{r) = S d ^ 2 f Pm {r)A{r)r d - 2 dr, (B8) 
J o 

which is exactly Ai(r) as defined in Eq. (j56j) . So, one has the identity M(r) = M.{r). 
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